Mixing effects for the structural relaxation in binary hard-sphere liquids 
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We report extensive molecular-dynamics simulation results for binary mixtures of hard spheres 
for different size disparities and different mixing percentages, for packing fractions up to 0.605 and 
over a characteristic time interval spanning up to five orders in magnitude. We explore the changes 
in the evolution of glassy dynamics due to mixing and discover two opposite scenarios: for large 
size disparity, increasing the mixing percentage of small particles leads to a speed-up of long-time 
dynamics, while small disparity leads to a slowing down. These results agree with predictions based 
on the mode coupling theory for ideal glass transitions. 
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This paper deals with structural relaxation in liquids, 
i.e. the dynamical phenomena which are precursors of 
the glass transition. These processes are characterized 
by rather complex dependencies of correlation functions 
on time and parameters like the density and the tem- 
perature. The introduction of several new experimental 
techniques and great progress in molecular dynamics sim- 
ulation studies has unfolded a wealth of facts about the 
evolution of structural relaxation upon cooling or com- 
pressing glass- forming liquids Q, @ • 

To establish an understanding of structural relaxation 
and the glass transition, one should focus on systems as 
simple as possible. Atomic one-component systems with 
conventional interaction potentials cannot be used since 
they crystallize before structural relaxation is fully de- 
veloped. Beginning with the work of Hansen and collab- 
orators |3| , binary mixtures of simple particles were used 
for molecular dynamics studies of the glass transition. A 
binary Lennard- Jones system has been used extensively 
in recent years to analyze structural relaxation Q . 

In previous work, mixing was introduced merely as a 
means of suppressing crystallization [|| . In the following 
we analyze systematically the effect of mixing on struc- 
tural relaxation to identify the influence of composition 
changes and variation of the particle size disparity on 
the glassy dynamics. We focus on binary hard-sphere 
mixtures (HSM), i.e. particles interacting via an hard- 
core potential, which we study via extensive molecular 
dynamics simulations. 

A specific motivation for our studies comes from the 
light-scattering data by Williams and van Megen for a 
HSM of colloidal particles [6j- The system was prepared 
to approximate a binary HSM for the ratio 6 = 0.60 of 
the particle diameters. Three mixing effects have been 
reported. If the percentage of the smaller particles in- 
creases from 10% to 20% of the relative packing fraction, 
then (i) the time scale for the final decay of the density 
correlators decreases; (ii) the plateau value for intermedi- 
ate times increases; (iii) the initial part of the structural 
relaxation slows down. The first effects means that mix- 



ing has stabilized the liquid as if the smaller particles 
provide some lubrification. This effect has some analog 
in the plasticization observed in dense polymeric liquids 
due to mixing with polymers of shorter lengths. However, 
the effects (ii) and (iii) indicate a stiffening of the dynam- 
ics upon mixing. These effects have never been reported 
for conventional systems, and one might wonder whether 
these are structural relaxation phenomena rather than 
colloid-specific features caused by, e.g., hydrodynamic in- 
teractions or polydispersity. Since our simulation studies 
are done implementing a Newtonian microscopic dynam- 
ics, this question will be answered by our results. 

Another motivation of our work is provided by recent 
predictions for binary HSM based on mode-coupling the- 
ory (MCT) calculations. This theory, which allows 
first principle evaluations of density-fluctuation correla- 
tors within a regime where structural relaxation dom- 
inates the dynamics, explains the light-scattering data 
for slightly polydisperse hard-sphere colloidal suspen- 
sions [8J. Extending MCT to binary HSM, the three 
mentioned mixing effects have been identified as struc- 
tural relaxation properties . The light scattering data 
for the 5 = 0.60 mixture could be described quantita- 
tively to a certain extent [Toj| . However, surprisingly, the 
theory predicts two different scenarios. The speed-up of 
the dynamics reported in Ref. [6j was found only for suf- 
ficiently large size disparity, say a size ratio S < 0.65. 
For 6 > 0.8, the opposite effect was predicted, i.e., mix- 
ing slows down the dynamics and the ideal-glass critical 
packing fractions decreases. The data reported in this 
letter confirms these astonishing predictions. 

We perform standard constant-energy molecular dy- 
namics |ll| simulation for a binary mixture of 1237 hard- 
sphere particles. The two species (A and B) have diam- 
eters c?a and de respectively, with d\ > de- The masses 
of the two species are taken as equal, thus all particles 
have the same thermal velocity, denoted as v t h- Units of 
length and time are chosen such that dx = 1 and vth = 1 . 
We use the size ratio S — d^/d^ as a control parameter 
specifying the size disparity. To model a large size dispar- 
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ity, we choose 8 = 0.60. To model systems of small size 
disparity we study the value 8 = 0.83. For each of the 
two 8 values, we study two different values of the relative 
packing fraction of the small species, x = tp-gli^A + Vb); 
where ip a — (n /6)g a d^. For the 5 = 0.60 system, we 
have studied x = 0.10 and x — 0.20, corresponding to a 
fraction of B particles equal to 34% and 54%. For the 
8 = 0.83 system, we have studied the cases x = 0.276 
and x — 0.37, corresponding to having 40% and 50% of 
B particles. For each 8 and x values, we study several 
values of the total packing fraction ip — pa + ^b, covering 
a region where dynamics slows down by 4 decades. Dur- 
ing the long simulation runs, we check for crystallization 
by monitoring the time evolution of the pressure and by 
visual inspection of the configurations. We also evalu- 
ate the structure factor to make sure that no crystalline 
peaks have developed. 

The ip dependence of the diffusion coefficient, evalu- 
ated from the long time limit of the mean square displace- 
ment 0| > is shown in Fig. ^ Note that the variation of 
the diffusivities extends over more than four decades. 

It is important to separate mixing effects already ob- 
served in normal liquid states from mixing effects which 
are peculiar of the structural relaxation dynamics. In 
hard-spheres, mixing effects for the normal-liquid dy- 
namics are qualitatively explained within Enskog's ki- 
netic theory. The diffusion constant D of the hard-sphere 
system (HSS) is expressed in terms of Enskog's collision 
rate v 1 D E — 3v^ h /(2i/). This rate modifies Boltzmann's 
collision rate for dilute gases by the contact value <?aa 
of the pair-distribution function: v = A^/Trv^ h gAgAAd A 
[l2j |. Enskog's formula is readily generalized to mixtures 
with equal constituents' masses, giving = 3v^/ (2v a ) , 
with v a — 4y / 7T«| h Vg Qpg a p$ a B- Here g a p denotes the 
pair-distribution function of the colliding pair a and (3 
for the distance at contact, d a p = (d a + dp)/2. Results 
for are included in Fig. [I] as lines for p < 0.39, using 
values from Percus-Yevick theory for g a p. The relevant 
normal-liquid mixing effects for the HSM, as given by the 
Enskog theory, can be summarized in the following: (i) 
Da decreases and Db increases upon decreasing 8; (ii) 
both Da and Db decrease upon increasing x; (iii) all Da 
values are smaller than all the Db values. A detailed 
discussion of normal mixing effects for an approximation 
to HSM has been published in Ref. |l3j . 

In agreement with the pioneer work of Alder et al. for 
the HSS the data for D a are described correctly on 
a 30% accuracy level by the Enskog theory for ip < 0.40. 
Upon increasing <p beyond 0.40, the cage effect starts 
dominating and leads to a suppression of D and D a below 
the corresponding Enskog values. Still we find that up 
to the large value <p = 0.48 there arc only normal-liquid 
mixing effects since the log D Q -versus-<^ curves do not 
exhibit crossings. 

Increasing p further, the curves in Fig.^get a stronger 
bending, i.e., the cage-effect enhances. If p reaches 
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FIG. 1: Diffusivities D a , a — A, B, for binary hard-sphere 
mixtures with size ratios 8 — de/rfA = 0.60 (open symbols) 
and 8 = 0.83 (filled symbols); packing contributions of the 
smaller spheres x = p-e.jp as in the legend, (a) Results for 
the two different 8 with a — A (circles) and a = B (squares) . 
The solid lines for tp < 0.39 are the Enskog results for Db(8 = 
0.6), D B (8 = 0.83), D A (8 = 0.83), and D A (8 = 0.6), from 
top to bottom. Dotted lines demonstrate power-law fits with 
exponents 7^(5, x), see text for details. The inset shows D* = 
fly 1 " to demonstrate the extrapolation to zero diffusivity. 
(b) Variation of Da upon changes of x for both 8 values; 
circles (triangles) refer to the higher (lower) x value studied. 



0.56, the diffusivities are about one order of magnitude 
smaller than in a conventional liquid. Within the interval 
0.48 < <p < 0.56, there occur crossings of the curves, i.e. 
the mixing properties change qualitatively. The cross- 
ing observed in Fig. ^ demonstrates a first remarkable 
feature of structural relaxation, viz., differently from En- 
skog predictions, diffusion of B particles for 8 = 0.83 
becomes smaller than Da for the 8 — 0.60 case. For 
if > 0.56, a new pattern evolves, showing the mixing 
effects for the fully developed structural relaxation. To 
support the association of the slowing down of the dy- 
namics to the approach towards a glass transition, the 
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FIG. 2: Normalized density-correlation functions for binary 
hard-sphere mixtures at a total packing fraction ip — 0.600 for 
particles with diameter ratio 5 = 0.60. The full and dashed 
lines present mixtures with relative packing contributions of 
the smaller particles B of x = 0.20 and x — 0.10, respectively. 
Upper and lower panel exhibit the results for wave vectors 
qd,A = 5 and qd,A = 10, respectively. The results for the 
particles B have been shifted down by 0.3 relative to the ones 
for the particles A. 



diffusivities have been fitted by a power-law function, 
D a = r a (<y,x)[^(5,x)-^]T«( 5 ' x ). As shown in Fig. □ the 
power-law lit can account well for the diffusivity changes 
over more than two orders of magnitude. For the four 
mixtures studied, we find that the critical packing frac- 
tion of species A deviates from that of species B by not 
more than 0.0008, so that we can confidently conclude, in 
agreement with MCT predictions, that each mixture ex- 
hibits a single critical packing fraction ip c (S,x). In the 
asymptotic limit ip — > ip c (S,x), the MCT exponent 7 
should be the same for Da and Db- This is not the 
case for our fit results, as also found in previous studies 
[lif . possibly because of preasymptotic correction effects. 

The presence of crossings in the diffusivity curves on 
entering the region where structural relaxation becomes 
relevant clearly shows that dynamics for 6 — 0.83 is sig- 
nificantly slower than the one for 5 = 0.60. Moreover — 



differently from the regime of normal-liquid dynamics 
where the diffusivities D a (8, x) exhibit the same trend 
with changes of x for 6 — 0.83 and S — 0.60 — there is a 
qualitatively different x-dependence of the long-time dy- 
namics for small and large size disparity, as recently pre- 
dicted by MCT 0. As shown in Fig. Eh while at S = 0.83 
dynamics becomes slower on increasing x, the opposite 
behavior is observed at S = 0.60. This second remark- 
able feature — the plasticization phenomenon alluded to 
before and exhibited by the 6 = 0.60 simulation data — is 
clearly shown in Fig^D. Consistently with these find- 
ings, at 6 = 0.60, the critical packing fraction increases 
on going from x = 0.10 {(p c {x = 0.10) = 0.6139 ± 0.0004) 
to x = 0.20 (f c {x = 0.20) = 0.6169 ± 0.0004). The crit- 
ical packing fractions for the system with S = 0.83 show 
the opposite trend: they decrease on increasing x from 
x = 0.276 (tp c (x = 0.276) = 0.5881 ± 0.0004) to x = 0.37 
{ip c {x = 0.37) = 0.5877± 0.0004). This third remarkable 
finding means that mixing stabilizes the glass for the sys- 
tem with small size disparity and destabilizes it for large 
size disparities. Different from the plasticization effect 
discovered by Williams and van Megen for S = 0.60, for 
6 = 0.83 an increase of the concentration of the smaller 
minority particles leads to a slowing down of the long- 
time density-fluctuation dynamics. The different scenar- 
ios for the long-time relaxation scales are demonstrated 
also in Figures |21 and [31 They show the x dependence of 
the density auto-correlation functions (f> a (q,t) for a = A 
and B at two different wave vectors, below (qd,A — 5) and 
above (qd,A = 10) the first peak of the structure factor 
(qd,A = 7), respectively. 

Within MCT, the ideal glass states are characterized 
by arrest of the density fluctuations. Within the liquid 
state, for ip below but close to f c , the cj> a (q, i)-versus- 
log t curves exhibit a plateau at (q) , the so-called crit- 
ical Debye- Waller factor. The curves are the flatter the 
smaller f c — ip, and the length of the plateau increases. 
This plateau, an outstanding feature of structural relax- 
ation, is clearly exhibited in Figs. [3 and |3J The whole 
structural-relaxation process consists of two steps. The 
first step deals with the relaxation towards the plateau 
fail)- The second step is the long-time process dealing 
with the decay of <f> a (q,t) from the plateau to zero. The 
three mixing features discussed above concerned the time 
scale of the second relaxation step. 

The results for (f>A(q,t) shown in Figs. [5] and for 
qd,A = 5 exhibit a fourth remarkable mixing effect: upon 
increasing x, the plateaus of the correlators increase. Ac- 
companied with this is a flattening of the </>a(9, i)-versus- 
logt curve in the plateau region. An indication of the 
same effect is also shown by the <^>b (q,t). These findings 
are in accord with the two mixing anomalies discovered 
in the experimental study of the colloidal S — 0.60 system 
0. Our results show that they are structural-relaxation 
phenomena that occur also in systems with Newtonian 
dynamics, and that they occur for small as well as for 
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FIG. 3: As in Fig. |2| but for a packing fraction tp = 0.582, 
size ratio 8 = 0.83, and relative packing contributions of small 
particles x = 0.37 (full lines) and x = 0.276 (dashed lines). 

large size disparity of the particles. We also note that 
on increasing x at 8 = 0.6, the increase of the height of 
the plateau value, associated with the speed-up of the 
dynamics, forces the relaxation curves for the same tp 
but for different x to cross at the beginning of the sec- 
ond relaxation step. For 5 — 0.83, an increase of x at 
fixed tp generates a slower dynamics and hence there is 
no crossing of the curves in the time interval for the sec- 
ond relaxation step. 

To conclude, we have studied structural relaxation of 
four binary hard-sphere mixtures via molecular dynamics 
simulations over dynamic ranges extending over about 
five orders of magnitude. Surprising mixing effects for 
the slow dynamics have been identified. In particular, 
we have shown that increasing the mixing percentage 
of the smaller minority particles can lead to a speeding 
up as well as to a slowing down of the long-time de- 
cay processes, depending on whether the size disparity is 
large or small, respectively. There is also an increase of 
the plateau of the density-autocorrelation functions for 
small and intermediate wave vectors reflecting a stiffen- 
ing of the nearly arrested glass structure. These findings, 
which pose a challenge to theories of the glass transition, 



show in particular, that the description of a glass-forming 
mixture by an effective one-component liquid cannot be 
possible for all properties of interest. The reported re- 
sults also confirm the conclusions arrived at in a light- 
scattering study of a quasi-binary colloidal suspension 

, and predictions within mode-coupling theory . Fi- 
nally, our work suggests that mode-coupling theory can 
contribute to an understanding of qualitative trends for 
the microscopic details in the dynamics of glass-forming 
liquids. 
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